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y— i ' Abstract 

o : 

\£) , We study inclusive charged-hadron production in collisions of quasireal pho- 

tons at next-to-leading order (NLO) in the QCD-improved parton model, using 
r S^ • fragmentation functions recently extracted from PEP and LEP1 data of e + e~ an- 

nihilation. We consistently superimpose the direct (DD), single-resolved (DR), and 
double-resolved (RR) 77 channels. We consider photon spectra generated by elec- 
tromagnetic bremsstrahlung and/or beamstrahlung off colliding e + and e~ beams 
as well as those which result from backscattering of laser light off such beams. First, 
we revisit existing single-tag data taken by TASSO at PETRA and by MARK II 
at PEP (with e + e~ energy y/S ~ 30 GeV) and confront them with our NLO cal- 
culations imposing the respective experimental cuts. We also make comparisons 
with the neutral-kaon to charged-hadron ratio measured by MARK II. Then, we 
present NLO predictions for LEP2, a next-generation e + e~ linear collider (NLC) 
in the TESLA design with y/S = 500 GeV, and a Compton collider obtained by 
converting a 500-GeV NLC. We analyze transverse- momentum and rapidity spectra 
with regard to the scale dependence, the interplay of the DD, DR, and RR com- 
ponents, the sensitivity to the gluon density inside the resolved photon, and the 
influence of gluon fragmentation. It turns out that the inclusive measurement of 
small-pT hadrons at a Compton collider would greatly constrain the gluon density 
of the photon and the gluon fragmentation function. 
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1 Introduction 



The inclusive production of single hadrons in fixed-target and colliding-beam experiments 
has been one of the most important testing grounds for the QCD-improved parton model. 
In contrast to the collective observation of jets of hadrons, the detection of single mesons 
and baryons also allows one to study fragmentation, i.e., the mechanism of how partons 
(quarks, gluons, and photons) turn into hadronic matter. In the framework of the QCD- 
improved parton model, the cross section of inclusive single-hadron production is described 
as a convolution of the parton-parton scattering cross sections with the parton density 
functions (PDF's) of the initial-state particles and the fragmentation functions (FF's), 
which characterize the transition of the partons that come out of the hard scattering 
to the hadrons that finally hit the detector. The factorization theorem |IJ ensures that 
the PDF's and FF's are universal and that only the partonic cross sections change when 
different processes are considered. While the partonic cross sections may be perturbatively 
calculated from the QCD Lagrangian, this is not yet possible for the PDF's and the FF's 
of hadrons with masses smaller than or comparable to the asymptotic scale parameter, A, 
and one has to determine them by fitting experimental data. The PDF's of protons and 
photons, which are needed to describe, e.g., pp, ep, and 77 reactions, are already highly 
constrained by measurements of deep-inelastic ep and e^ processes. 

The most direct way to obtain information on the FF's of hadrons is to analyze their 
energy spectrum measured in e + e~ annihilation, where the theoretical predictions are 
not obscured by additional nonperturbative input, e.g., in the form of PDF's for the 
incoming particles. After the extraction of leading-order (LO) FF's for pions and kaons 
from low-energy data of e + e~ annihilation and muon-nucleon deep-inelastic scattering in 
the late 70's and early 80's 0, there had long been no progress in this field. Some time 
ago, new next-to-leading-order (NLO) FF sets for charged and neutral pions and kaons 
as well as for eta mesons were fitted to e + e~ data generated with well-established Monte 
Carlo (MC) programs, which are fine-tuned so as to describe well a broad selection of 
experimental data ||. An alternative approach is to directly fit to experimental data. 
We took this avenue for charged pions and kaons ||, |5| and neutral kaons || by analyzing 
SL AC-PEP and recent CERN-LEP1 data. The assumption that the s, c, and b (d, c, and 
b) quarks fragment into charged pions (kaons) in the same way, which we had to make in 
Ref. |§], could be discarded in Ref. thanks to the advent of accurate data from the 
ALEPH Collaboration at LEP1 0, in which the fragmentation of gluons, b quarks, and 
light flavours into charged hadrons was distinguished. In the meantime, these data have 
become available in published form ||, where the fragmentation of quarks into charged 
hadrons is also reported for an enriched c-quark sample. This latter information was not 
at our disposal at the time of our analysis ||. However, we verified that our c-quark 
FF's for charged hadrons are approximately in agreement with this new measurement. 
In order to probe the scaling violation predicted by the Altarelli-Parisi equations ||, 
we made comparisons with other e + e~ data collected at different centre-of-mass (CM) 
energies. To test the validity of the factorization theorem for fragmentation, we also 
confronted measurements of inclusive single-hadron production by HI [10| and ZEUS 



2 



nj at the DESY ep collider HERA and by UAl |TJ| at the CERN SppS collider with 



the corresponding theoretical predictions based on our FF's H || [13]] f] In all cases, we 
found good agreement. Existing data on inclusive single-hadron production in collisions of 
quasireal photons offer yet another opportunity to quantitatively check the factorization 
theorem, which we wish to seize in the following. 

The purpose of this work is to use our new FF's for charged pions and kaons || to make 
predictions for inclusive single-charged-particle production in 77 collisions, which can soon 
be confronted with first experimental data from LEP2. With the exception of inclusive 7r° 
production fl4| , there exist no NLO predictions for inclusive particle production at LEP2 
in the literature. We also consider 77 physics at a next-generation e + e~ linear collider 
(NLC) with CM energy v^S* = 500 GeV. Furthermore, we compare our NLO analysis 
with 77 data at low CM energies measured some time ago with the TASSO detector at 
DESY PETRA pf and with the MARK II detector at SLAC PEP @. These data are 



at rather low px, where our NLO formalism is not expected to yield reliable results, so 
that a meaningful test of our FF's may not be possible. 

As is well known, three mechanisms contribute to the production of quarks and gluons 
in 77 collisions: (i) In the production through direct photons (DD), the two photons 
directly couple to the quark lines in the hard-scattering amplitudes. At least to LO, no 
spectator particles travel along the photon axes, (ii) If one of the photons splits into a flux 
of quarks (and gluons), one of these quarks may directly interact with the second photon. 
The remaining quarks (and gluons) build up a spectator jet in the direction of the split 
photon (single- resolved (DR) 7 contribution). The 77 cross section of this mechanism 
depends on the PDF's of the photon. (Hi) If both photons split into quarks and gluons 
(double- resolved (RR) 7 contribution), two spectator jets appear. Since the PDF's of 
the photon take large values — at small x through the gluon part and at large x through 
the so-called point-like quark part — the DR and RR contributions are numerically of the 
same order as the DD one. We perform a consistent NLO analysis, i.e., we include the 
DD, DR, and RR hard-scattering cross sections, the photon PDF's, and the FF's at NLO. 

The outline of this work is as follows. In Section 2, we shortly describe the formalism 
which we use to calculate the DD, DR, and RR contributions up to NLO. Section 3 
contains our numerical results for PETRA, PEP, LEP2, and NLC energies. Here, we 
also compare with the TASSO and MARK II data. Our conclusions are summarized in 
Section 4. 



2 Formalism 

When we speak of 77 collisions with almost real photons, we have in mind the e + e~- 
collision process where the positrons and electrons act as sources of nearly massless, 
collinear photons, which collide with each other to produce a debris of hadrons in the 
final state. The energy spectrum for the produced photons is well described by the 

^he NLO results shown in Rcf. Q were taken from Ref. |||, where the FF's of Ref. |§ were 
employed. 
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equivalent-photon approximation (EPA). In this approximation, the longitudinal- and 
transverse-momentum components decouple. The transverse momentum is integrated 
out with certain constraints, specified by the experimental situation. The resulting spec- 
tra for the photons emitted from the positrons and electrons are thus functions of the 
longitudinal- momentum fractions, x\ and x 2 , respectively. These distributions depend 
on the experimental setup. In our analysis, we consider past and future experiments at 
the colliders PETRA, PEP, LEP2, and NLC (TESLA design or laser spectrum). The 
respective forms of the EPA will be specified later when we present our numerical results. 
Our notation is as follows. We consider the reaction 

e + (pi) + e-(p 2 ) -> e+(pi) + e~(j/ 2 ) + h(p h ) + X, (1) 

where h is the observed hadron and X includes all unobserved hadrons. The four- 
momentum assignments are indicated in the parentheses. In the EPA, the inclusive cross 
section of process ([[]) is related to that of the corresponding 77 reaction, 



through 



lM)+l2(pl)^KPh)+X, (2) 



_ d 3 a(e + e — ► e + e h + X) } , , , 

E h — = J dxx I dx 2 F 7l/e +(xi)F 72/e -(x 2 ) 



d 3 a( 7ll2 -+h + X) 

x Eh rfv ' (3) 

where Xi = Ej/Ei (i = 1,2) and F li / e ±(xi) stand for the photon-spectrum functions to 
be specified later on. The lower bounds of integration, x ^ 2 j tire fixed by kinematics in 
terms of the transverse momentum p? and CM rapidity y of h. In the applications to 
follow, the actual e + and e~ acceptances must be incorporated in the integration ranges 
of x\ and x 2 , respectively. 

In the Q CD-improved parton model, the cross section of reaction (§) is expressed as 
a convolution of the parton-parton scattering cross sections, where one or both of the 
incoming partons may be photons, with the scale-dependent PDF's and FF's, 



d 3 0-(7i7 2 ->• h + X) fdx h ^^ ( _ A/f2 



J x z h V 



+ dx a F ah (x a ,M^)p\ 

a •* 

+ E/ dx b F bh (x b} M*)p] 



rf 3 q(7i72 -> c + X) 
d 3 p c 

Q d 3 a(a^ 2 — >■ c + X) 



d 3 p c 

d 3 a( ll b^ c + X) 
d 3 p c 

4~ ^ ] J dx a F a j^{x a , M 7 2 ) j dx b F bh (x b ,M*) 

a,b 

r .d 3 a(ab — > c + X) 
XPc *f c 



(4) 
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Here, the parton indices a, b, c run over the gluon and Np flavours of quarks and anti- 
quarks, k a = x a p\, k b = XbP2, and k c = Ph/xh are the parton momenta, F a / 7 (x a , M 2 ) is the 
PDF of parton a inside the photon, and Dh/ C {xh, M 2 ) is the FF of parton c into hadron 
h. The factorization scales, M 7 and Mh, will be specified later on. The first term on the 
right-hand side of Eq. (|j) is the direct-photon contribution (DD), the second and third 
ones are the once-resolved contributions (DR), and the fourth one is the twice-resolved 
contribution (RR). In LO, the DD contribution is of 0(a 2 ), where a is the electromagnetic 
coupling constant, i.e., it is of the same order as the the respective hard-scattering cross 
sections. The LO hard-scattering cross sections in the DR and RR components are of 
0(aas) and 0(a|), respectively, so that, at first sight, the LO DR and RR contributions 
might be considered to be of higher orders in the strong coupling constant, as- However, 
they are formally and numerically of the same order as the LO direct part, since the quark 
PDF's of the photon are enhanced by the factor a /as and the gluon PDF is enhanced at 
small x |T7|. In fact, the factor a /as from the photon PDF's properly adjusts the orders 



of the DR and RR contributions. 

It is convenient to define PDF's for finding a particular parton a, with momentum 
fraction x, inside the positron or electron, F a / e ±(x, M 2 ). In the EPA, this function is 
given by the convolution of the respective photon PDF, F a / y (x, M 2 ), with the photon- 
spectrum function, F y / e ±(x), introduced in Eq. (||), viz. 



■ ll ,±yx,M 2 ) = { d ^F ah (-,M 2 )F^ /e± (y). (5) 



If the photon directly participates in the hard scattering, the photon PDF in Eq. (|5|) 
must be replaced by the delta function 8(1 — x/y). Using definition fl5|), we can combine 
Eqs. (Q) and (||) in one formula. The inclusive cross section to NLO is then written as 

„ d 3 a(e + e~ — > e + e~h + X) r dxh „ , „ ,os „ , „ ,os ^ , , ^ 

E h ^ -= L = £ / d Xl dx 2 -^F a/e ±{ Xl ,M 2 )F b/e ±{x 2 ,M 2 )D h/c (x h ,M 2 ) 

d Ph „ i. ~J x h 



1 

X 

7TS 



-r — ( s , v ;v )5(l-w) 

v dv 



(6) 



where it is understood that the QCD-corrected hard-scattering cross sections in the DD, 
RR, and two DR channels are properly multiplied with the respective e + and e" PDF's 
and summed over. As usual, v = 1 + t/s and w = —u/(s + t), where s = (k a + kb) 2 , 
t = {k a — k c ) 2 , and u = (kj, — k c ) 2 are the Mandelstam variables at the parton level. 
They are related to the external Mandelstam variables, S = (pi +P2) 2 , T = (pi — P2) 2 , 
and U = (j>2 — Ph) 2 , by s = X\X 2 S, t = x{F/xh, and u = x 2 U/xh- The functions 
Kk a k b ^k c (s,v,w; fi 2 , M 2 , Ml) contain the NLO corrections to the hard-scattering cross 
sections, and fi is the QCD renormalization scale. 
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The Kk a k b ^k c functions for the DD channel have been derived by Aurenche et al. [ITS 



and have recently been confirmed by Gordon | l4fl . The DR Kk a k b ^k c functions have been 
calculated in Ref. [19] and reevaluated in Ref. [lj], where explicit expressions may be 
found. The RR Kk a k b ^k c functions have been obtained in Ref. ||20|| . They have recently 
been applied to predict single-charged-particle production in low-Q 2 ep scattering |2TJ . 
Direct and resolved photoproduction in ep scattering correspond to the DR and RR 
components in 77 scattering, and we may convert our previous analysis by replacing the 
proton PDF's with the photon PDF's. For the evaluation of the DD contribution, we 
employ a computer code created by Gordon in connection with Ref. [14] . 



For consistency, one needs both the photon PDF's and FF's in NLO. As photon PDF's 
we use the NLO set by Gluck, Reya, and Vogt (GRV) f22|, which we translate into the 
MS factorization scheme. In their analysis, the M" 2 evolution starts at a rather low value, 
Mq = 0.3 GeV 2 . Therefore, this set is also applicable at rather small pt- However, one 
should keep in mind that the predictions may not be reliable in the small-py region because 
as is large and nonperturbative effects may dominate. As FF's for charged particles, i.e., 
the sum of charged pions and kaons, we employ our recently constructed NLO set ||, 
which has been extracted in the MS factorization scheme. This scheme has also been 
employed for the derivation of the NLO kernels Kk a k b ^k c - 

In the calculation of the Kk a k b ^k c functions pertinent to the DD and DR channels, one 
encounters collinear singularities associated with the splitting of the incoming photons into 
qq pairs. These initial-state singularities are absorbed into the bare PDF's appearing in 
the resolved-photon contributions, i.e., the DR and RR components, which renders these 
functions M 7 dependent. A similar M 7 dependence, but with opposite sign, shows up in 
the Kk a k b ^k c functions of the DD and DR processes. In this way, the DD, DR, and RR 
processes become interrelated. The DR (RR) part is actually a NLO contribution to the 
DD (DR) part that cannot be treated fully perturbatively anymore. This nonperturbative 
part is then described by the photon PDF's. Up to higher-order terms in the photon 
PDF's, the M 7 dependence cancels in the combination of the DD and DR contributions 
on one hand, and the DR and RR contributions on the other hand. It is clear that the 
classification into DD, DR, and RR contributions becomes ambiguous if NLO corrections 
are included. NLO terms that have been attributed to the DD (DR) contribution may be 
shifted to the DR (RR) part. The cancellation of the dependences on M 7 and the choice 
of factorization scheme associated with the incoming-photon leg in the superposition of 
direct and resolved photoproduction in ep scattering has been demonstrated in Ref. |2l| . 
We expect that this cancellation also works for the 77 process (fj). The formal aspects of 



such cancellations in 77 reactions have been investigated by Gordon [14 



3 Numerical Results 

We are now in a position to present our numerical results for the cross section of inclusive 
charged-particle production in 77 collisions at NLO. Our plan is as follows. First, we 
consider existing data taken at low energies by TASSO |T5j and MARK II jnjj and confront 
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them with the corresponding NLO calculations. Then, we make predictions for LEP2 
and NLC, superimposing the energy spectra of the usual EPA bremsstrahlung off the 
incoming positrons and electrons and of beamstrahlung according to the TESLA design 
to be discussed later. As a further application, we consider 77 collisions with photon 
spectra as they result from Compton backscattering of laser light off the e + and e~ beams 
of the NLC. If this method works, it will be possible to generate high-luminosity beams 
of real photons carrying approximately 80% of the e + (e~) energy (see Ref. |2t| for the 
impact of various NLC collider options). 

We work in the MS renormalization and factorization scheme with Np = 5 active quark 

(5) 

flavours and Ai^ = 131 MeV, except for our TASSO and MARK II analyses, where we use 

Np = 4 and = 200 MeV. We identify the various scales and set /1 = M 7 = Mh = £pt, 
where £ is a dimensionless scale factor. Unless stated otherwise, we put £ = 1 and employ 
the NLO set of the GRV photon PDF's @ in the MS scheme. We take the NLO FF's for 
charged pions and kaons from our recent work |J. We sum over charged charged pions 
and kaons. Charged-baryon production is likely to be negligible; in e + e~ annihilation, 
some 10% of the charged tracks are due to protons and antiprotons ||], and one expects 
their fraction to be even smaller in 77 processes. In our MARK II analysis, we calculate 
the ratio of the Kg and charged-hadron cross sections. For this purpose, we employ the 
K s FF's from our very recent work ||. 

In the following three subsections, the results are ordered according to the various 
energy regions relevant for PETRA/PEP, LEP2, and NLC. 



3.1 Results for PETRA/PEP Energies 

Before presenting our predictions at high energies (LEP2 and NLC), where we expect 
to see experimental data at large pp in the near future — at least from LEP2 — , we wish 
to investigate how well the QCD-improved parton model can explain existing data at 
lower energies. Unfortunately, there are only few data for the inclusive cross section with 
information on the absolute normalization. Such data come from TASSO at PETRA 



|15| and MARK II at PEP (T6|. In the TASSO experiment, the effective CM energy was 
Vo = 33.1 GeV, which is the average for runs with beam energies between 13.7 and 
18.3 GeV JSJ]. The tagging requirements for the positrons and electrons were such that 
one of the two leptons (positron or electron) was tagged in the forward detector, which 
covered a narrow angular region between # min = 24 mrad and 9 max = 60 mrad relative 
to the beam direction, if it carried an energy of at least 4 GeV, whereas the other lepton 
(electron or positron) remained untagged jyj. In our theoretical evaluation, we take these 
tagging conditions into account by calculating the energy spectrum of the photons from 
the tagged leptons by means of the formula 



7 7 /e(x) ^ 



1 -hi ^) 2 ln QL £ + 2 2 



x ™ O 2 e \ O 2 




(7) 



where m e is the electron mass, Q 2 min = (1 — x)E 2 9 2 nin , and similarly for Q 2 max (E = vS/2 
is the beam energy). The tag-energy condition leads to a photon-energy cutoff at x max = 
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0.76. The no-tag photon spectrum is calculated from the well-known formula by Brodsky 



et al. [25|, written in a convenient form in Ref. [26]. Furthermore, we include an overall 
factor of 2 to account for the fact that both the positron or the electron can trigger the 
forward detector. The average Q 2 of the tagged photon was (Q 2 ) = 0.35 GeV 2 [|l5j. In 
our evaluation, we take Q 2 = in the hard-scattering cross sections and photon PDF's. 
This is justified since (Q 2 ) <C p\ for pr<^2 GeV, which we are primarily interested in. 
A further experimental constraint in the TASSO experiment was | cos 6*| < 0.84, where 



9 is the angle of charged tracks with respect to the beam axis flip] . This restricts the 
y range over which we integrate the doubly differential cross section d 2 ajdydp\ to be 
\y\ < 1.22. Taking these kinematical constraints into account, we calculate da /dp 2 -, for 
charged particles. 

The result is shown in Fig. |1] as a function of pt for the three different scale choices 
£ = 1/2, 1, 2 and compared with the TASSO data. As we can see, the experimental cross 
section is systematically larger than the theoretical one. In particular at larger px, where 
our approach is supposed to be valid, the measurement greatly exceeds our prediction, 
even for £ = 1/2. At small px, we do expect that the data overshoot our theoretical result, 
since in this region soft interactions give rise to additional particle production, which has 
not been subtracted from the data. This contribution is supposed to be significant for 
px^l-5 GeV and is usually described as in hadronic reactions induced by the vector- 
dominance-mechanism (VDM) component of the photon. We have no explanation for 
the disagreement with the data in the upper px range. A serious background for the 
production of high-p^ events from one-photon annihilation events, where either of the 
incoming positron or electron radiates a hard photon, was thoroughly investigated by 
TASSO |15j and was subtracted from the data shown in Fig. [1| It is interesting to study 
the relative importance of the DD, DR, and RR contributions. At px ^ 1 GeV, the RR 
component is dominant, whereas for larger px the bulk of the cross section comes from 
the DD channel. At px = (2, 3, 4) GeV, the DD component amounts to (51, 69, 75)%, 
respectively. This shows that, at these px values, the use of different photon PDF's would 
not change the cross section in any appreciable way and that the disagreement with the 
data is due to the DD component. 

Single-tag measurements similar to those by TASSO were performed with the MARK II 
detector at PEP [16| operating with y/S = 29 GeV. The virtualities of the photons 
emitted from the tagged lepton were required to lie between Q 2 min = 0.075 GeV 2 and 
Qmax = 1 GeV 2 and had the average value (Q 2 ) = 0.5 GeV 2 fllBI . The x range of these 
photons and the y range of the produced hadrons were only constrained by kinematics 
||16|| . The MARK II data are confronted with our theoretical predictions for £ = 1/2, 1, 



and 2 in Fig. 0. The theoretical curves are multiplied by an overall factor 1.2, to account 
for px smearing and resolution effects not corrected for in the experimental data [27j. We 
observe that the agreement between theory and experiment is considerably better than in 
the case of TASSO. At rather low px, where theory is not expected to be reliable because 
of the importance of the soft hadronic component, which was not subtracted in the data, 
there is almost perfect agreement, while the theoretical curves undershoot the data at 
large px- The dent at px ~ 3 GeV in the curve for £ = 1 is due to the onset of c-quark 
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fragmentation [0 . The disagreement between the MARK II charged-particle cross section 
and the theoretical predictions at larger px has previously been noticed by Aurenche et 
al. [^8] and by Gordon [14]. However, these authors used older FF's, partly in LO f25fl . 



Our results, which are based on FF's that are rigorously constrained by recent e e data 
0, suggest that this disagreement does not originate from the FF's. It also appears that, 



at large px, our results exceed those of Refs. [14, 28 



The MARK II Collaboration also measured the cross section of inclusive Kg produc- 
tion under kinematic constraints identical to those used for their charged-hadron sample 



i~6|] . Using our recently constructed set of Kg FF's ||, we calculate the ratio of K° 



s 



to charged hadron production as a function of px- The results are compared with the 
MARK II data in Fig. [|, showing satisfactory agreement for px ^ 1 GeV. The step in 
the theoretical curve at px ~ 3 GeV comes about because the c-quark threshold affects 
charged-hadron and Kg production somewhat differently. 



3.2 Results for LEP2 

Unfortunately, there are no LEP1 data available which could be directly compared with 
our theoretical results. ALEPH [29] and DELPHI |3(| reported on inclusive measure- 
ments of charged tracks, extending up to px ~ 4 GeV. However, these are not absolutely 
normalized, and the EPA constraints, which are indispensable for absolute predictions, 
are not specified. 

It is hoped that better data, which extend to even larger px, will become available 
very soon from LEP2. For our predictions, we assume v^S* = 175 GeV and describe the 
quasireal-photon spectrum in the EPA by the formula |31| 

2^2 



a \l + (l- x )\ E 2 9 2 max (l - xf + mix 



F l/e {x) = — { In 

+ 2(1 -x) 



2ir [ x mix 2 

mix 1 



E 2 6 2 n Jl — x) 2 + m 2 x 2 x 



We vary x = E 1 /E e over the full range allowed by kinematics and put the antitagging 
angle to 6 max = 30 mrad. 

In Fig. [|, we show the DD, DR, and RR contributions to the cross section da/dpx of 
inclusive charged-hadron production and their sum as a function of px- As in the case of 
MARK II, da/dpx is obtained from d 2 a /dydpx by integrating over the full y range. For 
Px^5 GeV, the cross section is dominated by the DD component, which is then bigger 
than the sum of the DR and RR contributions. In general, the DR and RR components 
yield only a small fraction of the total sum, except for very small px [pr ^ 3 GeV) , where 
the RR part dominates. The y spectrum for px = 10 GeV is plotted in Fig. |^, where 
the DD, DR, and RR contributions are again displayed along with their sum. Due to the 
symmetric experimental setup, the curves are symmetric in y. For px = 10 GeV, the DD 
contribution dominates for all values of y. It is of interest to know whether the resolved 
contributions can be used to obtain information about the gluon PDF of the photon. 
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Since for large pr the DR and RR contributions are small, this should be feasible only 
for small px- To investigate this point, we repeat the computation of da/dpx putting 
Fg/ 1 {x) = and plot the ratio of the outcome to the full result in Fig. |6] (dashed line). We 
observe that this ratio is below 80% only if px ^ 5 GeV. This means that this cross section 
is not a very powerful discriminator to disentangle the gluon PDF of the photon. Almost 
the same pattern is observed if the gluon FF is switched off (dot-dashed line). The cross 
section da/dpx is sensitive to D g j h ±(x) only for small px- This may be understood by 
observing that the gluon FF falls off much more rapidly with increasing x than the quark 
FF's. In Fig. ||, we also show the fraction of da/dpx originating from the fragmentation 
into K . It is rather large, ranging between 34% and 40% for pt between 5 and 30 GeV. 

It is generally believed that the scale dependence of a parton-model calculation to 
a given order indicates the size of unknown higher-order corrections and thus may be 
used to estimate the theoretical uncertainty of the prediction. This leads us to study in 
Fig. |7| the £ dependence of da/dpr and its DD, DR, and RR parts for p T = 10 GeV. This 
scale dependence is determined mostly by the DD component, which dominates the cross 
section at px = 10 GeV. In the range 1/2 < £ < 2, the total sum varies by ±8% relative 
to its value at £ = 1. The discontinuity at £ = 0.95 stems from the 6-quark threshold. 



3.3 Results for NLC 



Next, we consider the predictions for the NLC with = 500 GeV and TESLA design. 
At the NLC, photons are produced not only by bremsstrahlung, but also via synchrotron 
radiation emitted by one of the colliding bunches in the field of the opposing bunch . 
This phenomenon has been termed beamstrahlung. The details of the beamstrahlung 
spectrum depend crucially on the design and operation mode of the NLC. In our study, 
we select the TESLA design, where the unwanted effects of beamstrahlung are reduced to 
a tolerable level. We coherently superimpose the EPA and beamstrahlung spectra. The 
EPA spectrum is computed from Eq. (|8]) with 9 max = 175 mrad and the beamstrahlung 
spectrum from the expression given in Ref. fl32"| |, with parameters T e ff = 0.039 and 
cr 2 = 0.5 mm |33]. This combined spectrum has recently been used also in a study of 
heavy-quark production in 77 collisions pl |. The pt distribution, integrated over the 
full y range, is shown in Fig. ^| together with its DD, DR, and RR components. Apart 
from an overall enhancement due to the increased EPA logarithm and the additional 
beamstrahlung-induced contribution, the pt spectrum looks very similar to the LEP2 
case in Fig. |[ The corresponding y spectra for px = 10 GeV are shown in Fig. |9|. Due 
to the admixture of beamstrahlung, their shapes differ somewhat from those of the pure 
EPA case in Fig. j|. 

The highest possible photon energies with large enough luminosity may be achieved 
by converting the NLC into a 77 collider via backscattering of high-energetic laser light 
off the e H 



and e beams [p5| . The resulting photon spectrum is given by |35] 
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where 

G M = fr-*-^W + «)4 + !-2uiKF (10) 



K K 2 



This spectrum extends up to x max — k/(1 + k). In our calculation we choose k = 4.83, 
so that x m ax = 0.83. If the experimental setup was arranged so that k > 4.83, e + e~ pairs 
would be produced in the collisions of the primary laser photons and the high-energetic 



backscattered photons. In Fig. [10], we present the DD, DR, and RR contributions to the 
Pt spectrum da/dpr together with their sum. Compared to the previous cases, the relative 
magnitudes of the DD, DR, and RR components have completely changed. At small px, 
at about 5 GeV say, the RR component is by far the largest, whereas the DD component 
is negligible. At large px, around 20 GeV, the relation is quite different; the DD, DR, 
and RR contributions all have the same order of magnitude. A similar situation was 
encountered in the analogous study of heavy-quark production fl3~4f . Consequently, in the 



lower pt range, the cross section da/dpx of inclusive particle production with a Compton 
collider will serve as a powerful tool to extract information on the photon PDF's. The 
corresponding y spectra for px = 10 GeV are plotted in Fig. O. As far as the relative 
importance of the various components is concerned, we recognize a pattern similar to 



Fig. [10], which refers to the integrated rapidity spectrum, at px = 10 GeV. On the other 
hand, the y spectra are quite different in shape and relative magnitude from those found 
in the TESLA case. This may be understood by observing that the laser-photon spectrum 
is peaked at the upper edge. The RR component is now dominant and rather flat; the DR 
component is somewhat smaller and exhibits a characteristic two-hump structure with a 
significant depletion in the central region; the DD component is greatly suppressed, but 
exhibits a shape similar to the TESLA case. As a consequence, there is a pronounced 
plateau in the y spectrum of the total sum. Inspired by the observations made in Fig. |Tl], 
we quantitatively study in Fig. |12| the influence of the gluon PDF of the photon. The 
dashed line represents the ratio of dajdpx with the gluon PDF switched off to the full 
result as a function of px- We see that, at p T = (5, 15, 25) GeV, (75, 42, 23)% of the cross 
section stems from the gluon PDF. The NLO photon PDF set by Gordon and Storrow 
(GS) P^] differs from the GRV set |22|] mainly in the gluon PDF. The ratio of da/dpx 



evaluated with the GS set to the GRV result is visualized in Fig. [12] by the solid line, 
which is a smooth interpolation. Obviously, the GS and GRV predictions appreciably 
differ at small px, by up to a factor of 2. Thus, the NLC operated as a 77 collider would 
provide an excellent laboratory to pin down the poorly known gluon PDF of the photon. 
In Fig. [12], we also assess the sensitivity of such an experiment to the gluon FF by showing 
the fraction of da/dpx due to the quark FF's (dot-dashed line). At px = (5, 15, 25) GeV, 
(45, 21, 13)% of the cross section is related to gluon fragmentation. Consequently, the 77 
NLC would also allow one to probe the gluon FF, which is only moderately constrained 
by data of inclusive hadron production via e + e~ annihilation. 
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4 Summary and Conclusions 



We studied the inclusive production of charged hadrons in collisions of almost real photons 
at NLO in the framework of the QCD-improved parton model. This approach is concep- 
tually very different from the one based on MC event generators, which is frequently 
followed by experimentalists to interpret their data. In these MC programs, the various 
QCD processes are simulated using LO matrix elements in connection with certain model 
assumptions concerning the formation of hadronic final states. Although these MC pack- 
ages often lead to satisfactory descriptions of the data, from the theoretical point of view, 
their drawback is that this happens at the expense of introducing a number of ad-hoc 
fine-tuning parameters, which do not originate in the QCD Lagrangian. Furthermore, 
in the MC approach, it seems impossible to implement the factorization of final-state 
collinear singularities, which impedes a consistent extension to NLO. On the contrary, in 
the QCD-improved parton model, such singularities are absorbed into the bare (infinite) 
FF's so as to render them renormalized (finite) in a way quite similar to the procedure 
for the PDF's at the incoming legs. In our analysis, we used FF's extracted with NLO 
precision from PEP and recent LEP1 data on inclusive single-hadron production in e + e~ 
annihilation H |6|]. In fact, owing to the factorization theorem PJ, the FF's are process 
independent as long as one considers unbiased single-hadron event samples. By the same 
token, our formalism is less flexible than the MC approach, since we are unable to describe 
more complicated final states, with more than one detected hadron and acceptance cuts 
on observables different from our integration variables. 

There are various ways to realize 77 processes in experiments. The conventional 
method is to exploit the QED bremsstrahlung which is radiated during e + e~ reactions 
from the colliding beams. High-quality data of this type were collected by TASSO |1| at 
PETRA and MARK II [[U| at PEP. Unfortunately, the bulk of these data are accumulated 
at rather small px, where the NLO formalism is not expected to yield reliable results, 
and the few data points at larger px values carry considerable error bars due to limited 
statistics. Detailed comparison revealed that our consistent NLO calculation with up-to- 
date FF's was able to reduce somewhat the mismatch between theory and experiment 
that had previously been observed by other authors [n, However, especially in the 
case of TASSO, the degree of discrepancy in the upper px range still gives some reason 
for concern. 

While at LEP1 the 77 processes were overwhelmed beneath the tremendous back- 
ground due to e + e~ annihilation on the Z resonance, the situation will be much more 
favourable at LEP2. We presented NLO predictions for the pt and y spectra to be mea- 
sured in the oncoming LEP2 phase. We verified that our NLO predictions are rather 
stable under scale variations; at px = 10 GeV, the cross section only fluctuates by ±8% 
if £ is varied between 1/2 and 2. Unfortunately, except at rather low p T , the study of 
77 processes at LEP2 will not much deepen our understanding of the nature of the gluon 
inside the photon and of its role within the fragmentation process. 

Photon-photon physics at the NLC will greatly benefit from beamstrahlung as an 
additional source of quasireal photons and, of course, from the increased EPA logarithm. 
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In the case of the 500-GeV NLC with TESLA architecture, the cross sections at pt = (5, 
15, 25) GeV exceed the corresponding LEP2 values by factors of (19, 26, 37), respectively, 
while the shapes of the pt and y spectra closely resemble their LEP2 counterparts. 

The advent of a Compton collider, generated by backscattering of laser light off the e + 
and e~ beams of the NLC, would represent an important step forward in the history of 
77 experimentation, and might considerably improve our knowledge of both the gluon FF 
and the gluon PDF of the photon. By converting the NLC, the cross sections of inclusive 
charged-hadron production at pt = (5, 15, 25) GeV would be increased by factors of 
(35, 52, 88), respectively. At the same time, the relative importance of the DD channel 
would be dramatically reduced, so that the sensitivity to the photon PDF's would be 
correspondingly amplified. 
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FIGURE CAPTIONS 



Figure 1: Cross section da/dpip of inclusive charged-hadron production in single-tagged 
77 collisions. The TASSO data JTS] (\/S = 33.1 GeV on average) are compared with our 
corresponding NLO calculations for scale choices £ = 1/2, 1, and 2. 

Figure 2: Cross section da/dp? of inclusive charged-hadron production in single-tagged 
77 collisions. The MARK II data [I6| (y/S = 29 GeV) are compared with our corre- 



sponding NLO calculations for scale choices £ = 1/2, 1, and 2. 

Figure 3: The ratio of inclusive Kg production to inclusive charged-hadron production 
as measured by MARK II [16] is compared with our NLO calculation for scale choice 
£ = 1. The rise near px — 3 GeV is related to the c-quark threshold. 

Figure 4: NLO cross section da/dpx of inclusive charged-hadron production in double- 
tagged 77 collisions at LEP2 (y/S = 175 GeV). The DD, DR, and RR components are 
also shown. 

Figure 5: Same as in Fig. | for the y dependence of d 2 a /dydpr at pt = 10 GeV. 

Figure 6: Influence of the gluon PDF of the photon, the gluon FF, and the charged-kaon 
final states on the £ = 1 result of Fig. Shown are the fractions that remain if the gluon 
is switched off in the photon PDF's (dashed line) or in the FF's (dot-dashed line) as well 
as the fraction due to charged-kaon production (solid line). 

Figure 7: £ dependence of the results shown in Fig. [| at pt = 10 GeV, in the interval 
0.3 < £ < 4. 

Figure 8: Same as in Fig. § for the NLC with TESLA design (y/S = 500 GeV). 

Figure 9: Same as in Fig. g for the NLC with TESLA design (y/S = 500 GeV). 

Figure 10: Same as in Fig. [| for the NLC in the Compton-backscattering mode (yfS = 
500 GeV). 

Figure 11: Same as in Fig. [5] for the NLC in the Compton-backscattering mode (y/S = 
500 GeV). 

Figure 12: Influence of the gluon PDF of the photon and the gluon FF on the total 
result of Fig. [TT]. Shown are the fractions that remain if the gluon is switched off in the 
photon PDF's (dashed line) or in the FF's (dot-dashed line) as well as the ratio of the 
calculation with the GS photon PDF's to that with the GRV set (solid line). The solid 
line is a smooth interpolation. 
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